Directed abelian algebras and their applications to stochastic models 
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To each directed acyclic graph (this includes some D-dimensional lattices) one can associate some 
abelian algebras that we call directed abelian algebras (DA A). On each site of the graph one attaches 
a generator of the algebra. These algebras depend on several parameters and are semisimple. Using 
any DAA one can define a family of Hamiltonians which give the continuous time evolution of a 
stochastic process. The calculation of the spectra and ground state wavefunctions (stationary states 
probability distributions) is an easy algebraic exercise. If one considers D-dimensional lattices and 
choose Hamiltonians linear in the generators, in the finite-size scaling the Hamiltonian spectrum is 
gapless with a critical dynamic exponent z = D. One possible application of the DAA is to sandpile 
models. In the paper we present this application considering one and two dimensional lattices. In 
the one dimensional case, when the DAA conserves the number of particles, the avalanches belong 
to the random walker universality class (critical exponent a T — 3/2). We study the local density of 
particles inside large avalanches showing a depletion of particles at the source of the avalanche and 
an enrichment at its end. In two dimensions we did extensive Monte-Carlo simulations and found 
a T = 1.782 ± 0.005. 
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I. INTRODUCTION 

It was Dhar [l|, [|J who pointed out that certain abelian 
algebras play an important role in understanding the 
properties of sandpile models. The time evolution of 
these systems is given, in the continuum limit, by a non- 
hermitian Hamiltonian which acts nonlocally in the con- 
figuration space. 

The expression of the Hamiltonian contains local gen- 
erators of an abelian algebra acting in the vector space 
of the regular representation of the abelian algebra. In 
some cases this vector space can be identified with the re- 
currence states (configurations) of the sandpile models. 
These algebras have a special structure describing top- 
plings from a site to its neighbors. The net result being 
a highly nonlocal interaction. 

The abelian algebra appears in two different ways. On 
the one hand it determines the Hamiltonian which drives 
the system to the stationary state. On the other hand, 
if one considers the action of one of the generators on 
the stationary state of the system, the algebra describes 
through topplings a " short" (instantaneous) time behav- 
ior of an avalanche triggered by the local action of the 
generator. One can attach a virtual time (this is the 
"short" time) to the evolution of the avalanche. This is 
possible due to the abelian nature of the algebras. In 
some cases the real and virtual times can get combined 
in a unique time (see [|| for an example), we are not 
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considering this case here. 

In the present paper we are going to present a large 
class of abelian algebras that we call directed abelian al- 
gebras (DAA). They are defined as follows. One takes 
a directed acyclic graph [4| (these graphs are much used 
in mathematics and computer sciences Q). They are 
formed by oriented links (edges) to which we attach ar- 
rows. There are no oriented closed loops. To each site 
(vertex) of that graph one attaches a generator of DAA. 
There are as many relations among the generators as sites 
in the graph. Examples of directed acyclic graphs are 
some lattices in any number of dimensions. The direc- 
tion in which the sandpiles propagate can be identified 
with the virtual time. The evolution of the avalanches 
can be seen as a growth process in one space-dimension 
less @. 

The DAA are semisimple (all representations are de- 
composable into irreducible representations). They also 
have a special property: they have only one-dimensional 
irreducible representations and the number of nonequiv- 
alent representations coincides with the dimension of the 
regular representation. Since the Hamiltonians act onto 
the regular representations, they can be diagonalized in 
a trivial way. This stays true even if the Hamiltonians 
describe the evolution of systems in any number of di- 
mensions. At least in some cases, all the eigenf unctions 
can be easily obtained. 

The stationary state, which is of product measure form 
(no correlations!), can also be obtained in any number D 
of dimensions. The thermodynamic limit of the spectra 
is also very simple. 

All systems described by directed algebras are critical 
with the dynamic critical exponent having the value z = 
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D. This does not imply that the avalanches have to be 
critical since they evolve in virtual time. They might 
have exponential tails. 

The paper is organized as follows. In Section 2 we de- 
fine a quadratic DAA on a one-dimensional lattice. It is a 
two-state model in which one has a particle or a vacancy 
on each site. We give the Hamiltonian which gives the 
time evolution of this system and show how to compute 
exactly its spectrum. The ground-state wavefunction is 
of product form. The average density is uniform in spite 
of having an open system. 

The Hamiltonian spectrum and eigenfunctions are 
given in Appendix A for the physically relevant case in 
which the number of particles is conserved. 

In Section 3 we study the avalanches in the model de- 
fined in Section 2. We show that they belong to the 
"random walker" universality class The probabil- 
ity to find an avalanche of virtual time T corresponds 
to the probability to have a first passage at time T of 
the random walker. The connection between thepresent 
paper and the totally asymmetric Oslo model @ [l(| is 
presented in Appendix B. We also look at the average 
density of particles left in the system when the avalanche 
stopped at time T. The density of particles is not uni- 
form anymore. The physical process being directed, the 
region near the beginning of the avalanche is depleted 
and the region at the end of the avalanche is enriched in 
particles. For long avalanches, the particle density dis- 
tribution functions become singular at both ends of the 
avalanches. It is easy to understand the physical origin 
of the two singularities. 

Higher order one-dimensional DAA are defined in Sec- 
tion 4. They correspond to N-state models. All the prop- 
erties of the 2-state models are recovered and no new 
qualitative features appear. 

The construction of a DAA on an arbitrary directed 
acyclic graph is presented in Section 5. As applications 
we consider DAA on a one-dimensional lattice, in which 
a site i is related to two successive sites i + 1 and i + 2, 
and on two-dimensional lattices. In two dimensions, one 
obtains avalanches belonging to the universality class de- 
scribed in pjl, [TH, 03] . We have done large scale Monte- 
Carlo simulations to measure the exponent a T (usually 
denoted by r f ). The result doesn't agree with the ex- 
pected value a T — 7/4. 

Our conclusions are presented in Section 6. 



II. A QUADRATIC DIRECTED ABELIAN 
ALGEBRA AND ITS APPLICATION TO A 
STOCHASTIC PROCESS. 

In this section we first show how to define the Hamil- 
tonian describing a stochastic process by using an asso- 
ciative algebra. We then consider a simple application 
of quadratic directed abelian algebra defined on a one- 
dimensional lattice. Associated to this algebra we may 
define a family of stochastic processes, sharing the same 



stationary state. 

Let us now consider an associative algebra with gen- 
erators A(i) (i = 1,2,..., to). Taking products of the 
generators one gets the words W(r). The algebra is de- 
fined by giving some relations between the words. If we 
can choose the linearly independent words {W(r);r — 
1, . . . , nw} such that any product of them verify the re- 
lation: 

nw nw 

W(r)W( S )=Y / P r g SW ( ( l)> P r q S >°> X>«' S = 1 > 

q=l q=l 

then the Hamiltonian 

H = J2<u)(l-W(u)), (2) 

u=l 

acting in the vector space defined by the basis of inde- 
pendent words {W(s)}, is an intensity matrix provided 
that the coefficients c(it) are nonnegative [la ]. H is act- 
ing from the left on the words W(s) of the vector space: 
W(u)W(s) (W{u) is one on the terms in © and W(s) 
belongs to the vector space). The action of the nw in- 
dependent words {W(w)} on the vector space defined by 
the same words gives the regular representation of the 
algebra. 

If the algebra A contains a left ideal [l6| defined by 
the words /, the Hamiltonian ((2]) acting on this ideal 
gives again a stochastic process. If A has several ideals 
I\ C I2 C • • • C /„, H has a block triangular form. The 
ground-state wavefunction (eigenvalue equal to zero) is a 
linear combination of the states (words) which define the 
vector space I\ . These states correspond to the recurrent 
states. The words belonging to the ideals I2, ■ ■ ■ , I n but 
not belonging to I\ do not appear in the stationary state. 

One can obtain Hamiltonians describing stochastic 
processes using associative algebras even if the positiv- 
ity conditions ([1]) are not fulfilled. This can be achieved 
by not working in the regular representation of the alge- 
bra but in other representations. The following simple 
example illustrates the procedure. 

Consider the quadratic algebra with only one generator 

a: 

a 2 = (1 — fi)a + fj,e, e 2 = e, ea — ae (3) 

and 

H = 1 — a, (4) 

fi is a parameter. If < n < 1, the condition ([1]) is 
satisfied. H is an intensity matrix (order e and a): 

»- (s -;)■ « 5) 

where e is the 2x2 unit matrix. Let us now perform the 
similarity transformation: 

H' = SHS-\C) (6) 
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where 



S = 



1 £ 
£ 1 



and £ = (1 - a)/{l-0). H' is: 



if' 



o 

—a 



-0 



(J) 



The algebra ((3]) stays unchanged with /j = at+ft—l. The 
matrix ([7]) is the most general 2x2 intensity matrix ( 
< a, (3 < 1 ). Notice that the physical process depends 
on two parameters instead of one (which appears in the 
algebra ((3])) and that the range of fj, has changed: — 1 < 
fi < 1. 

The probabilities P q (t) to find the words W(q) at the 
time t can be obtained from the master equation 



where 



-P(t) = -HP(t), 



P(t) = J2P q (t)W(q). 



(8) 



(9) 



We now apply this formalism to a simple example of a 
quadratic DAA (for applications to non-abelian algebras, 
see [Hj]). Consider a one-dimensional lattice with L sites. 
We attach to each site i (i = 1,2, ... ,L) a generator a,i. 
A quadratic DAA is defined by the relations: 



a \ = Pl a l+l +P2di+1 



for i = 1, 2, . . . , L — 1, 



H + (1 - /i)a£, 



[aj,aj] = 0, i,i = 1,2,...,L, 
where the p m and g n are probabilities and 

H=Pl+P2+P3, 1-^ = 91+92- 



(10) 



(11) 



(12) 



(13) 



Notice that the words with a given value i are related to 
words with values j > i, therefore we have a preferred 
direction on the lattice. The relation (flTj) was obtained 
taking a^+i = 1 in equation (| 10[) . 

The 2 L monomials 1, a\, . . . , Ol, aa.a.2, . . . , ol-iOl, ■ ■ ■ , 
aia2 ■ ■ ■ cll are the independent words of the DAA (they 
define the regular representations of the DAA). One can 
check, using (|10 p - (fT3")) , that the relations ([T|) are veri- 
fied. The DAA ([Tt) )) -([T2" |) can therefore be used to define 
a family of stochastic processes whose Hamiltonians are 
obtained by arbitrary choices of c(u) (u = 1, . . . , 2 L ) in 
©. If fi = 0, the algebra (fTt ty -lfTg] ) has an ideal 



h = aia 2 



• a L , 



(14) 



which gives in the regular representation an absorbing 
state in the stochastic process defined by ©. 

One can map the vector space of the regular represen- 
tation onto a lattice with L sites. On each site i one has 
a particle if appears in the monomial or a vacancy if 
di doesn't appear in the monomial. To the word " 1" one 
associates the empty lattice. There are 2 L configurations 
obtained in this way. 

The a-i (i = 1, 2, . . . , L) acting onto the regular repre- 
sentation are not in the diagonal form. They can how- 
ever be diagonalized simultaneously and one obtains in 
this way 2 L one-dimensional irreducible representations 
of the algebra. The eigenvalues of at denoted by Xi can be 
obtained recursively using the distinct eigenvalues a^+i 

Of Oi+l 



XL 
™2 



= 1,M) 

Xi(qix i+ i +q 2 ) - {pix 2 i+l +p 2 x i+ i +p 3 ) = 0, 



where i = L — 1, L — 2, . . . , 1. 

Notice that one has to solve only quadratic equations 
and that the <Zj have a common eigenvalue xi = 1. One 
can easily find the corresponding eigcnfunction in the 
regular representation: 



a 4 $ = $, 
The expression of $ is: 

$ = 



1,2,..., L. 



(15) 



(16) 



In the expression of $ enters only the parameter jj, instead 
of the four parameters appearing in (|10p . The proof is 
straightforward. We have: 



aL- 



\i + a L [i + a L 



1 



1 



(17) 



Where we have used (TlTj) . We now use (ITU1) and (fTTjl to 

obtain: 



O-L-l- 



( x + a L -i (i + cil _ m + a L-i M + a L 
l + /i l+/i 1 + n 1 + yu 

and repeat the procedure to get, due to (fl"2j). 

= i = L — 2, . .., 1. 



(18) 



(19) 



The eigenstate $, due to (fTBj). gives the stationary state 
for the arbitrary stochastic models defined by @, i. e. 
HQ = 0. The average density in the stationary state (fTB]) 
is equal to 1/(1 + fi) and varies between 1/2 and 1. 

Other equivalent representations of the algebra (flQ|) 
can be used to describe stochastic processes. In Appendix 
B we give two examples. One of them is the totally asym- 
metric Oslo model @, [To| ■ 

We now consider the simple stochastic process given 
by the Hamiltonian: 



H 



L 



(1 - (n)/L. 



(20) 
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The physical interpretation of the stochastic process de- 
scribed by (|20p can be more easily understood using the 
configuration space and discrete times. At the time t one 
takes a configuration with a certain distribution of par- 
ticles on the L sites (not more than one per site). With 
a probability 1/L a particle is added at the site i. If the 
site is empty, the particle stays at i. If one has a parti- 
cle at i already, one obtains several configurations with a 
probability given by the toppling rules given by the DAA 

Notice that the stationary state is of product form and 
therefore there are no correlations in the stationary state. 

To compute the spectrum of H, we can take all the 
in their diagonal representation. Since H acts in the reg- 
ular representation, all the irreducible one-dimensional 
representations of the DAA have to appear. The number 
of irreps coincides with the dimension of the vector space 
given by the regular representation. Therefore, each rep- 
resentation appears only once and the calculation of the 
spectrum is trivial (see Appendix A). It is interesting to 
compute the energy gap of H. The ground-state energy 
is zero. The first excited state is obtained taking <2j = 1 
(i = L, L— 1, 2). It follows from (JTUJ) that a\ can take 
two values: 1 and — fj,. One obtains two one-dimensional 
representations. The first one is the ground state with 
energy zero and the second one, corresponding to the first 
excited state with energy E\ = (l+fi)/L (we call energies 
the eigenvalues of H). The dynamic critical exponent z 
follows: 

lim LEi = 1+ fi, (21) 

L — >oo 

therefore z = 1. This result is interesting for several 
reasons. The value z = 1 is also obtained when the sys- 
tem is conformal invariant [Til ] . This is certainly not the 
case here since there are no correlations in the stationary 
state. The finite-size scaling limit of the spectra can be 
changed (through the change of the parameters in (fT0|) ) 
almost at will and are not described by representations of 
the Virasoro algebra as expected in a conformal invariant 
theory. The result is not surprising since in the configu- 
ration space, each at acts through topplings in a highly 
nonlocal way and there is no reason to expect conformal 
invariance. It is nice to have an example which illustrates 
this phenomenon in full generality. 

There is another interesting consequence of (f2"Tj) . Since 
the energy gap in the thermodynamic limit vanishes, the 
system is critical. Does this imply that the avalanches, 
to be discussed later, are described always by power laws 
as expected by SOC? The answer is no. As one can see 
from (|10p. unless all the coefficients P2,Pz and qi van- 
ish, the number of particles is not conserved, particles 
are lost. One expects therefore an exponential decay for 
the probability of having an avalanche of a certain size 
[]~7L [l8j . For avalanches which are critical one defines a 
dynamic critical exponent z a [l^ | related to the virtual 
time in which the avalanche occurs. It is not clear to us 
why the two exponents z and z a should coincide although 
in many examples they do. 



The spectrum and eigenfunctions of H are given in 
Appendix A for the simpler case in which one conserves 
the number of particles in (|10[) . The knowledge of the 
spectra allows us to get information of time dependent 
quantities. If one starts, for example, with an empty 
lattice, one can look at the time dependence of the av- 
erage density which, in the large L limit, should reach 
algebraically its value in the stationary state. 

There is a purely algebraic way to compute the density 
for finite L. Let us assume that at t = we start with an 
empty lattice. This corresponds to the word "1" in the 
algebra. In the basis of monomials (see ©), one has 

P(t) = e- Ht • 1 = e- Ht . (22) 

After expanding the exponential in (f2"2")l we use the al- 
gebra (fT0|) and expresses P{t) in terms of monomials. We 
take a,i — y (for all i's) and obtain the function P(t,y). 
The density is given by the expression: 

Pit) - ^ ■ < 23 > 

Let us observe that we can take the parameters in 
(fT0|) and (TTTj) z-dependent, the ground-state wavefunc- 
tion maintaining the product form: 

* = ]It£p (24) 
i=i 

where the definition of [n is obvious. The dynamic criti- 
cal exponent z stays unchanged. 

III. A QUADRATIC DIRECTED ABELIAN 
ALGEBRA AND AVALANCHES. 

In the last section we have discussed the time evo- 
lution of a system described by H given by ([2D)) . In 
order to study its spectrum and obtain the ground-state 
wavefunction we have used the quadratic DAA (TlT))) - (ri'5)) . 
From these results, in this section, we need only the ex- 
pression and use the DAA to study avalanches. To 
simplify the calculation we consider the case in which the 
number of particles is conserved. In this case, relations 
(TIT)]) take the form: 

a? = /xa? +1 + (1 - fjL)aia,i + i. (25) 

The physical meaning of (|25p is simple: if one has two 
particles on the site i, with a probability // they reach 
the site i + 1 and with a probability (1 — /i) one particle 
moves to the site i + 1 and one stays at the site i. Mul- 
tiplying (f2"5| by a power of a.; we can find what are the 
probabilities that n particles on site i move in block on 
the site i + 1 or only n — 1 particles move and one stays 
at the site i. Obviously the number of particles is not 
conserved at the end of the system. 
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To examine avalanches in our system, we take the sta- 
tionary state given by (fT!)|) and consider only the config- 
urations in which the first site is occupied (they have a 
probability 1/(1 + fi)). We add one particle on the first 
site and analyze the consequences. The overall effect of 
adding the particle is trivial to derive: 



L L 

2 TT I- 1 ' a i / / \ TT I 1 ' a i 



I' 



(26) 



where we have used (fT5)l and p^|) . This tells us that the 
changes take place only at the first site. With probabil- 
ity fi it becomes empty and with probability (1 — (a) is 
occupied by a particle. 

We are interested in more detailed properties of this 
process. First we would like to know what is the proba- 
bility p(k) to have the first k sites affected by the action 
of a\ and the remaining L — k sites untouched (this de- 
fines an avalanche of size k, starting at the first site). The 
physics inside an avalanche is easy to understand looking 
at (|25|) . Particles close to the origin of the avalanche are 
pushed towards its end. We will be interested therefore 
to obtain the non-homogeneous density of particles inside 
an avalanche. Let us first obtain p(k). 

We use (|2"5j) and obtain: 



i=1 



\xa\ 



1 + /it 



1 



/' 



i=3 



/' 



(27) 



One can see that we have gotten an avalanche given by 
ai 02 of size 2 (the first two sites occupied) with proba- 
bility jtt(l — A t )/(1 + A 1 )- If we are not interested in the 
content of the avalanche, and only on its size, we can 
formally take a\ = 1 in (|27p to get: 



M + a % - A*(l - A*) a 2 + (1 — A* + A* 2 ) a 2 + / ia 2 



i=2 
L 



i=3 



/' 

(J, +JH 

A 4 ' 



1 + /J 



(28) 



If we want to proceed further, and compute the proba- 
bilities of larger avalanches we have to compute quantities 
like: 



1 + A* 



(29) 



-jj x a t a i+1 \, 



in which we have taken into account the fact that the 
algebra (|25|) conserves the number of particles. Multi- 
plying (|29[) by ai one obtains the recurrence relations: 



B, 



(n+1) _ ,,ri{n). n (n+l) _ R (n) 



C ( n+1} = liD^'; D^ 1 ' = Cg 1 ' + (1 - aW(30) 



AtCp; C{ ' = B K ' + (1 - Ai)Q 



(n) 



with 

cf = 1, D<°> = 0, Cf > = 1 - M , £>« = /x(31) 
From which we get: 

£>(") = /iC^"" 1 ) (32) 

and 

c (n+l) _ c (n) = _^ C {n) _ C [ n -V). (33) 

The solution of the recurrence relation (l33l) is 



(n) _ 1 - (-M) 



n+1 



1 + /X 



(34) 



In order to look for avalanches, in (|29p one takes a« = 1 
and using (f3T)|) and (j3"2")l , one obtains 



t Ai + a; 



1 



1 + A* 1 + M 

+ {c^+^c { r 2) K +1 +^ { r 1] ^]- (35) 

If A* 7^ 1 (for fx = 1 the process is deterministic and the 
two particles on site 1 move on site L where they leave 
the system) from ([54]) one can see that for large values of 

n, C[ becomes independent on n and equal to 1/(1 + /x). 
Substituting this result in (|35[) one obtains: 



„ A* + Oi 



1 



1 + A* (l + /i) ; 



[Kv + i+(i+/)«r + i+KVi]- (36) 



One could have obtained this last equation without 
solving the recurrence relation (f3"Tj) - (|33[) (we are going to 
use this method in the next section). Equation (|33[) has 

(n) 

an n-independent solution C\ = C. We introduce this 
solution in (|35|) and ask for the sum of the coefficients to 
be equal to 1. This gives C = 1/(1 + li). 

The whole virtual time evolution (r > 1), of the 
avalanches is given by the following expression 



L t L 



71=1 



j = T+l 



/' 



(37) 



where P n (j) is the probability that at virtual time r an 
avalanche, with n particles at site i = r, is taking place. 
These probabilities, from (|25|) and (|35|) . satisfy the recur- 
rence relations 

Pi{r) = P 2 (T-1)RV\ 

F»(t) = P 2 (r - 1)4 2) + P 3 (r - l)i? W , 

P„(r) = P„ +1 (r-l)i? ( ! l+1) +P„(r-l)4 n) 



for 2 < n < t, with 



(38) 
(39) 
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and 



tin — 



f(n-l) 



R 



(n) 



1 - f (n - 



n+1 



1 



.(40) 



We can visualize these last equations in the following 
way. We take two coordinates: the virtual time r (dis- 
critized values) as the horizontal axis and a coordinate 
n (discritized values) as a vertical axis. Equation (|37p 
describes a random walker which at time r is at the posi- 
tion n, and moves to the positions n — 1 or stays at 
n with probabilitoes which are n-dependent. The proba- 
bility of having an avalanche of size k is given by the first 
passage at virtual time r = T = k to return to its initial 
position x = 1, and is given by p(r) = Pi(r). This prob- 
ability is not simple to calculate due to the n-dependence 



of the random walk parameters R_ , R, 



(n) D (n) 



n 



and R^ in 



(|38|) . However the large size avalanches are dominated by 
random walks that stay mostly at large distances from 
the origin (n large). We take = + fi) in (l4"0"l) 
and obtain 



= R = fi/{l+iif, R ( n) = 1-2R = 



1 + fJ, 2 



(1 



(41) 

With probabilities /i/(l + fi) 2 , the walker moves to n + 1 
or n — 1 and with a probability (I + fi 2 ) / (1 + /i) 2 stays at 
n. In the continuum, one obtains the diffusion equation 
with a diffusion constant D = /i/(l + /i) 2 . Consequently, 
the probability of having large avalanches of size k cor- 
responds to the usual first passage time at t = T = k 



p(T)=P 1 (T) 



1 



VDT* 



(42) 



If one defines the exponent a T by the relation p(T) 
1/T CTt , where T is the duration of the avalanche, one 
obtains a T = 3/2. The approximation (|4"Tj) for the prob- 
abilities (|4"0|) is valid if n >> 1/| ln/x|. Since n is of the 
order of y/T, we conclude that one can see the avalanches 
described by the power law (|4"2"]l for 



T»T C = l/(ln/z) 2 



(43) 



In order to illustrate the approximation (|4Tj) we com- 
pare in Fig. [T] the results obtained by the exact recur- 
rence relations with fi — 1/4 ([38]) (black curve) with 
the one obtained by using the asymptotic relations (|4"Tj) 
(red curve). We only show the probabilities for relatively 
small avalanches (r < 50) where the curves show a more 
pronounced difference (T c = 1 in this case). We also show 
(blue dots), for comparison, the average results obtained 
from 2 x 10 4 samples in a Monte Carlo simulation. 

Having in mind applications to two-dimensional sys- 
tems (see Section 5) we are interested to know how large 



~i 1 i r 

Monle Carlo -2. 10 4 samples 
First passage (exact) 
First passage asymptotic 



b 0.04 - 



I 

20 
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FIG. 1: (Color online) Exact and Monte Carlo results for 
the probability distribution of avalanches of size T for the 
stochastic model defined by J20]) and (J25J) , with fi = 1/4. The 
dots are the results from Monte Carlo simulations. The black 
curve is the exact result, calculated algebraically. The red 
curve is obtained by considering the asymptotic behavior of 
the recurrence relations (I41|l in (|38[) . 



have to be the avalanches to allow us to get a good esti- 
mate for <t t . We consider the quantity: 



o r {T) = T 



p(T )-p(T+l) 
P(T) 



(44) 



which is an estimate for a T . In Fig. @ we show the val- 
ues of ov(T) for fj, = 0.99 (T c = 9900 in this case). We 
observe that <r T (T) first increases with T giving a maxi- 
mum value of 1.91 and then decreases. Notice that even 
taking very large avalanches one still gets a poor esti- 
mate: <7 T (10 6 ) = 1.5445. This estimate is larger than the 
expected value 1.5. Using extrapolation techniques (8- 
degree polynomial fitting) we can improve substantially 
the estimate and get <r T = 1.5006 ± 0.0008. We have to 
stress that for p(T) we have used exact results and not 
Monte Carlo simulations as in Section 5. 

If (i is small, the estimates of cr T (T) improve dramat- 
ically. In Fig. ^ we show o>(T) for \i = 0.1 (T c is 
negligible in this case). 

We observe that in contrast with the case fi = 0.99, the 
estimates increase with T and one can get a very good 
precision for relatively small values of T. A best fit to 
the data gives : 



p(T) = A/T 9 + B/T h , 



(45) 



where A = 1.0794 ± 0.0003,5 = -7.16 ± 0.05, g = 
1.50000 ± 0.00001, h = 2.448 ± 0.003. 

In this way one has obtained a concrete realization of 
the model of Maslov and Zhang 0] in which z a = 1. The 
same result was obtained in a different context. A gen- 
eralized totally asymmetric hopping processes on a ring 
(fixed density of particles) was considered More 
than one particle per site are allowed. The rules of the 
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2 ~ T= 15,326 




1.6 - 
1.5 - 

14 l , 1 , 1 , 1 , J 

0.0001 0.0002 0.0003 0.0004 

1/T 

FIG. 2: The estimate a T {T) defined in gU) for ^ = 0.99 and 
2400 < T < 10 6 . The estimate increases, with T, up tp T = 
15,456. The last value of the estimate is <t t (10 6 ) = 1.5445. 




1/T 

FIG. 3: The estimate a T (T) defined in (O for p = 0.1 and 
1 < T < 10 6 . The last value of the estimate is a T (10 6 ) = 
1.49999. 



hopping processes are fixed by the condition of integra- 
bility of the system (one can write the Bethe ansatz). 
Interestingly enough the rules coincide with those ob- 
tained from the DAA (|25|) . There is no relation (fTTjl on 
a ring. If the density of particles on the ring is equal to 
the average density 1/(1 + p) obtained in our model, the 
two models describe the same large avalanches. 

We would like to make the connection between our ap- 
proach and the work of Stapleton and Christensen on 
the totally asymmetric Oslo model (TAOM) and its gen- 
eralizations [10| | . where similar results were obtained. As 
shown in detail in Appendix B in the case of the two-state 
model, the TAOM corresponds to take another represen- 
tation (not the regular one) of the same algebra (f25|) . 
This representation depends on two parameters whereas 
the algebra depends only on one. The same picture 
stays valid for the multi-state one-dimensional models 



discussed in Section 4. We have directly proven the cor- 
respondence to the random walker process and proceeded 
by deriving the probability distribution of avalanches. In 
10] the authors have computed moments of the toppling 
probability distribution and came to the same conclusion. 

Up to now we were interested in the probability of 
having an avalanche of size T. We would like to know 
what is the physics of a given avalanche. It is useful to 
visualize the virtual time r either as a time coordinate 
or a space coordinate. 

Looking at the DAA (f23|) we suspect that particles 
present in the stationary state (average density < p >= 
1/(1 + p)) are pushed away from the origin of the 
avalanche towards its end. For a given avalanche of size 
T we consider the quantity : p(r,T)— < p >, where 
p(r, T) is the avergae density at virtual time r. In order 
to compute this quantity we could use ([29]) without tak- 
ing di = 1. This is a lengthy calculation. Instead we did 
Monte-Carlo simulations. In Fig. [4j for p = 1/2 we show 
(p(T,T)- < p >) multiplied by ^T/D as a function of 
t/T for three avalanche sizes. The function is antisym- 
metric. D is the diffusion constant (D = p/(l + p) 2 )- 
We notice not only a nice data collapse but a singular 
depletion of particles for small values of t/T and a sin- 
gular enrichment of particles for t/T close to one. We 
have done similar simulations for other values of p and 
the data come on top of those shown for p = 1/2. If one 
looks r as a space coordinate, one gets what one would 
expect to see in the aftermath of an avalanche. 

For small values of t/T one finds: 

c (D\ 1/2 
P-<P>-— A [-) (46) 

where c ~ 0.56. 



4 



2 



-4 





M.C. - T=250 N av . 


= 5.10 4 


1 1 












— M.C. - T=500 N >v> 


= 2.5 10 4 






M.C. - T=1000 N aja 


= 5 10 3 






Predicted T = 1000 








i , 














,1 


1,1, 



_J I I I I I l_ 

0.2 0.4 0.6 0.8 

T/T 



FIG. 4: (Color online) Depletion and enrichment of particles 
in avalanches of different sizes T. (T/ D) 1/2 (p(r,T)- < p >) 
as a function of t/T obtained from Monte-Carlo simulations 
for T = 250, 500 and 1000. The numbers of avalanches N ava 
observed in each case are given in the insert. The solid curve is 
the result of the calculation using (gSJ) and igSJ for T = 1000. 
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Theoretically one can obtain the scaling function which 
describes the data shown in Fig. 2] in the following 
way. We use the random walker approximation. Let 
P(x, r; T) = P(x, T — t;T) be the conditional probabil- 
ity to find a random walker at time < t < T at the 
position x if the first passage time is T. This is a Brow- 
nian excursion Q. The average position at time r is 



< x > T ,T= / xP(x,T;T)dx. 
Jo 

One can compute 

1 d 

P{T,T)- < p > = --- < X > T . T 



(47) 



(48) 



In order to do the calculation we have used the lattice ex- 
pression P(x, t; T) obtained in the simpler case in which 
the random walker which is at the position x at time r 
moves at x + 1 or x — 1 with probabilities 1/2 at r + 1. 
The random walker starts at x = at time r = and 
returns at x — at r = T (T even). One has [2l| : 



P(x,t;T) : 



x 2 T{{T/2- l)!) 2 



2 (T- 
l)!(T-r- 



-2)! 



(49) 



Note that in the continuum the diffusion constant is 
D = 1/2 in this case. The function P(x,t;T) satisfies 
various scaling laws. In particular if we take t = T/2, it 
scales as shown in Fig. [5] A consequence of this scaling 
behavior is that 



< X >T 



T i/2_ 



(50) 



For t = T/2, < x > gets its maximum value. One has 
checked a more general scaling property: 



< x > TlT - 



fDTf(r/T). 



Obviously D = 1/2 in the present case. 
For small values of r, one gets: 



< x >- 



(51) 



(52) 



where c = 0.56 . The expression (|52[) is what one expects 
for a random walker. The constant c appears from the 
constraint x > 0. If we use (gSJ) and 02$ we get (gBJl, 
therefore the singularity observed at small values of r has 
a simple explanation. 

We have computed < x > t ,t and used l|48p to obtain, 
in this way, the solid curve shown in Fig. 3J 




FIG. 5: (Color online) The conditional probability function 
P(x,T/2;T) multiplied by T 1/2 as a function of x/T 1/2 for 
T = 500, 1000 and 2000. 



As we are going to see all the results obtained for the 
quadratic case generalize in a trivial way. The order N 
DAA with conservation of the number of particles is de- 
fined as follows: 



N-l 



-.N 



= J2 a ka^ + - 1 k ,(i = l,...,L-l), (53) 



k=0 



N-l 



/ J otha L , 

k=0 



N-l 



^2 a k = 1, < a k < 1. 



(54) 



(55) 



fc=0 



One recovers the case N = 2, taking a = A 4 - The sta- 
tionary states $ of the family of Hamiltonians (see ©) 
associated to these algebras, i. e., 



a,$ = 1, 



have the expression: 



Em=0 



LI a j s^N-1 



(56) 



(57) 



where 



IV. HIGHER ORDER ONE-DIMENSIONAL 
DIRECTED ABELIAN ALGEBRAS. 



0m = y^, a k- 

k=0 



(58) 



The quadratic DAA discussed in the previous section 
can be generalized to higher orders. These algebras de- 
fine stochastic models with recurrence states (stable con- 
figurations) containing at most N — 1 particles per site. 



The average density 



(59) 



9 



varies between (N — l)/2 and N — 1. 

The algebra has A L one-dimensional representa- 
tions. As in section 2 we restrict ourselves to the simple 
Hamiltoman ([20]) . Its eigenspectrum can be obtained in 
the same way as for the N — 2 case. As compared with 
the case N = 2 one has to solve iV-th order algebraic 
equations instead of quadratic ones. The dynamic criti- 
cal exponent stays unchanged 2 = 1. 

We present the study of the avalanches for the case 
N = 3 only. The generalization to other values of A*" will 
become obvious. 

We take cto and oti as the parameters of the algebra. 
Using (|57|) we have 



a i + ( a o + otx)a,i + a 



2ao + ai 



(60) 



We consider only configurations in which on the first site 
we have two particles and add one more particle on this 
site. As in the case N = 2, we have to compute the 
quantity: 

aj + (op + ui)ai + a x 
ai l + 2a + «i = ( 1 + 2a o + a 1 ) 



(61) 



Multiplying (|6ip by Oj, one obtains the following recur- 
rence relations, valid for n > 2, 



We notice that (|65|) has n-independent solutions valid 
for large n: 



cf > = a 



(66) 



This implies (see (|63])) 



4 



(n) 



„ - a C, = a (a + ai)C, C^' ; = a^C, 



B[ n) = (a + ai)C, Cf J = (a + ai) a C, 



(n) 



(n) 



a (ao + ai)C. 



(67) 



If we are not interested in what happens in the 
avalanche on the site i, we take di = 1 in ([61]) . We 
substitute the coefficients in (IrTTI) by (j6l)|) - (f6"T|) and fix the 
value of C = (1 + 2qq + £*i) by imposing that the sum 
of the coefficients is equal to 1. We obtain 

a, — — = (1 + 2a + ai) 

1 + 2a + ai 

x [a a"+i 2 + (ao + «i)(l + aoW+Z 

+ (l + a2 + ( ao + a 1 ) 2 ))a™ +1 

+ (a + ai)(l + aoX+i 1 + a a a 1+i] ■ (68) 

This can be interpreted as a random walker which moves 
in a symmetric way one and two steps up and down. We 
are therefore back to the Maslov Zhang 0] universality 
class. 



DIRECTED ABELIAN ALGEBRAS IN 
VARIOUS GEOMETRIES. 







= -a (C 2 


-ct 


- 2) ) 






-cii(C 2 


-ct 








/ n( n ) 




- 2) ) 






1 7-1 (") 


-Df 




£,(«+!) 




— -«o(-C'2 


-Et 


- 2) ) 








- Et 





(62) 



A(n+1) 
A 




= ai c { 2 n) - 


f a Cf " 1} 


R (n+1) 




= ai D { 2 n) 




°0 


= a Q E { 2 n \D[ n+1) 


= ai E { 2 n) - 


f ao^ n_1 \63) 



Up to now we have considered DAA on a one- 
dimensional lattice. DAA can be generalized to more 
dimensions and different lattices. For simplicity, in this 
section, we confine ourselves to quadratic DAA which 
conserve the number of particles. 




V 

B 



with 



1 



E. 



(2) 



1 + 2a + a\ 
a 



£j(2) _ Qq + a\ 
' 2 ~~ 1 + 2a + ai ' 

, Cf } = = = 0.(64) 



1 + 2a + ol\ 
From we obtain 

-^2 — (,"o + aiJG 2 ,i/ 2 — a C 2 



(65) 



FIG. 6: The bulk site b is connected with the sites 1,2, ... ,M. 
Bulk sites are denoted by empty circles and boundary sites 
B by filled circles, respectively. Squares denote sites that can 
be of bulk or boundary type. 

We divide the sites of a graph into two sets: bulk sites 
and boundary sites. Take a bulk site b on the lattice and 
join it by outgoing arrows to the sites i (i = 1, 2, . . . , M) 
(see Fig. [B]) . To each of the M + 1 sites one associate a 
generator a^. A quadratic directed relation can be asso- 
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FIG. 7: Example of a directed acyclic graph, where the empty 
(filled) circles are bulk (boundary) sites. 



dated to the site b 

M 



<ij, 



k=l 

We denote 

N 



M M 

^(a k al + f3 k a b a k ) + ^ ^ Jk,ia k ai. (69) 

k=l l=k+l 



M M N 
fc=l fe=l/=fe+l k—l 



k, 



(70) 



where < /x < 1. To a boundary (sink) site B (no 
outgoing arrows) one associates a generator as satisfying 
the relation: 



a% = fi + (1 - fi)a B - 



(71) 



The sites i = 1,2, ...,M can be bulk or boundary 
sites to which we associate relations like (|69]l or (17T1) . 
respectively (see Fig. [7|) . There are as many relations as 
sites. The graph starts with bulk sites and ends with 
boundary sites. This is a directed acyclic graph 

We call a lattice with C sites on which one can define 
a DAA, a directed acyclic lattice. The DAA is defined 
by the £ relations and ([71]) • The algebra has 2 C 
one-dimensional irreducible representations and the reg- 
ular representation is 2 C dimensional. Examples of 2-d 
directed acyclic lattices are shown in Fig. [8] If we are 
interested in critical exponents we have to consider well 
defined sequences of larger and larger graphs. 

One can show that for any directed acyclic lattice one 
has: 



a/$ = $, I E £ 



where 



Jec 



[i + a j 



(72) 



(73) 



/ and J are all the sites of the directed acyclic lattice 
with C sites. The proof is by finite induction, similar to 
the one used for the one-dimensional case. One starts 
with the boundary generators and finishes with the bulk 
generators attached to sites with no incoming arrows. 



One can define a simple stochastic process on the di- 
rected acyclic lattice by the Hamiltonian 



H = 1 



C 



or, 



(74) 



and the C quadratic relations (|6"9"|) - (fTTj) . Obviously the 
stationary state of the system is given by (|T3"|) : 



H$ = 0. 



(75) 




FIG. 8: Examples of 2d directed acyclic lattices. In (a) the 
lattice is periodic in the horizontal direction, being defined on 
a cylinder. In (b) and (c) the lattice is open in the horizontal 
direction. 

In order to find the spectrum of H, one uses the diag- 
onal representations of the generators of a/. The eigen- 
values of aj are obtained using the one-dimensional rep- 
resentations of the DAA. One starts with the boundary 
generators (this is the last row in the lattices of Fig. [7]), 
solving the equations (|7Tj) . The solutions are introduced 
in the quadratic relations (|69|) related to the sites in the 
next layer (this is the row before last in Fig. [8]). The layer 
by layer procedure is repeated (we go against the arrows) 
up to when we have exhausted the lattice. The energy 
gap corresponding to the first excited state, is obtained 
taking for all the generators the eigenvalue 1 except for 
one generator belonging to the top of the lattice (there is 
only one in Fig. [5Jd), and many in Figs. [8^) and[8j;)) for 
which we take the eigenvalue — fx. This gives an energy 
gap: 



£?i = (l + n)/£, 



(76) 



and we conclude that the system is always critical and 
if the dimension of the directed acyclic lattice is D, the 
dynamic critical exponent is equal to z = D. 

Before discussing 2-dimensional directed acyclic lat- 
tices, we consider again a one-dimensional case in which 
the bulk sites in ([69]) are joined to two successive sites. 
The DAA is: 



a 2 + /3 2 a L _i + ja L 
/i + (1 - n)a L 



[3 2 aiai 



+2 



aia| 



(77) 
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We have shown that in general, we spare the reader 
the effort of having to follow a long proof, the avalanches 
are again given by a random walker. We would like to 
mention only an interesting observation. Let us take all 
the coefficients (l77|) zero except for 7. We get the deter- 
ministic algebra: 



af = a i+ ia i+2 (i = 0, 1, . . . , L) 



2 

°>L-1 



and the stationary state is (see ([75]) ) 



(78) 



(79) 



For simplicity, we take L + l sites (i = 0, 1, . . . , L) with 
L even and assume that we have two particles on the site 
i = which trigger the avalanche. We also group in the 
ground-state wavefunction (|79[) the sites in pairs. Using 
([75)1 and ([731) we nave 

L/2 

= aia 2 J|(l + a 2i -i){l + a 2 i)/4 



_ ,ai«2 a 2 a 3 a4 (1 + a 2 )a§a4 

~~ \ 5 A A / 



L/2 



n 

i=2 
L/2 

n 



(1 + a 2 j-i)(l + a 2 j) ^,1 03^4 
4 ^4 4 

i=2 



a|a4 s 



(1 + <22i-l)(l + 02i) 



(80) 



There are three terms in the last equation. The first 
one gives with a probability 1/4 an avalanche of size 2. 
The second term, with a probability 1/4, has the same 
structure as the same equation for a lattice with L — 2 
sites. Using the identity: 

2 (l + 02i_l)(l + 02i) 
a 2i-l a 2i ; 



(1 + a 2i _i)(l + a 2l )a 2i+ ia 2l+2 „ 2 
- 4 -a 2i _ia 2 . i+2 ,(81J 

we find that the last term in ([8"U]l describes 3 parti- 
cles which "fly" to the boundary and leave the system. 
This occurs with probability 1/2. Notice (see ([75])) that 

To sum up, if L is large, one gets with a probability 
l/4 fc an avalanche of size 2k. This implies that with a 
probability 1/3 one gets finite avalanches. With a proba- 
bility 2/3 one gets "avalanches" of the size of the system. 

Let us mention that in two dimensions a similar deter- 
ministic model gives avalanches with long tails (exponent 
o T = 3/2) P. 

We have studied the behavior of avalanches in the two- 
dimensional directed cyclic lattice shown in Fig. [SJd using 
the quadratic DAA algebra 



°>i,j = ot{na 2 i+hj + (1 - n)a^a i+1}J ) + 
(1 - a){na 2 i+1 A + (1 - n)a id a i+ltj ). 



The avalanches observed in the lattice in Fig. [8j) coin- 
cide with those observed in the cylinder in Fig. [5Ji if the 
perimeter of the cylinder is larger than the size of the 
avalanche. 

We did extensive Monte Carlo simulations for several 
values of the parameters in the DAA. We took only small 
values of /1 to make sure that we don't get crossover ef- 
fects like those described in Fig. O We have considered 
up to 1.2 x 10 7 avalanches of sizes up to T = 30000. It is 
important to note that the estimate is smaller than the 
expected value. 

The exponent a T was obtained in several ways. A pos- 
sible estimate is the value of ui — a T that makes 



dF(u,T) 



dT 



for 17, 500 < T < 30, 000, where 



F(uj,T) =T U 



= 0, 



p(t)dt. 



(83) 



(84) 



(82) 



In Fig. (O we show for a = 1/2 and /i = 1/2 the func- 
tion F(uj,T) for some values of 10. The figure indicate 
1.77 < a T < 1.78. The same calculation were done for 
a = 1/2, /i = 1/10 giving us 1.780 < cr T < 1.785, and also 
for the anisotropic lattice where a = 1/4, fi — 1/2 giving 
us 1.78 < o> < 1.79. Our results indicate that the model 
in two dimensions belong to a universality class where 
a T = 1.78 ± 0.01. We should keep in mind that in one 
dimension for [i = 0.1, and T = 30000, one has obtained 
an estimate <J T (T) = 1.49971 , quite close to and smaller 
than the correct value 3/2. Previous Monte-Carlo sim- 
ulations [IH had larger errors. If one takes our results 
at face value they contradict the theoretical predictions 
[HI, Q3] ■ We cannot exclude however that larger values 
of T are necessary to get the ultimate value of oy . 



VI. CONCLUSIONS 

Using directed acyclic lattices Q we have defined di- 
rected abelian algebras (DAA). The lattices can be in any 
dimensions. We attach a generator of the DAA to each 
site of the lattice. There are polynomial relations among 
the generators. A given DAA can depend on several pa- 
rameters. The remarkable property of these algebras is 
that they are in general semisimple. 

For each DAA one can define a family of Hamiltonians 
acting onto the regular representation of the DAA. These 
Hamiltonians give the time evolution of various stochas- 
tic processes. The vector space of the regular represen- 
tation can be mapped in the vector space space obtained 
taking N states on each site of the directed acyclic lat- 
tice (N is the degree of the polynomial relations). In this 
way one obtains iV-state Hamiltonians in an arbitrary 
number of dimensions. 

It is easy to compute the whole spectrum for the family 
of Hamiltonians related to a given DAA. All one has to 
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FIG. 9: (Color online) The function F(u, T) as defined in (JS4J) 
for some values of 10. This figure is for the stochastic model 
associated with (|82p with a — 1/2 and fj, — 1/10. In the inset 
we show the linear fit for the values of u> obtained from the 
datas 17, 500 < T < 30, 000. 

do is to find all the one-dimensional irreducible represen- 
tations of the algebra. This is a purely algebraic exercise. 
Once the spectrum of the generators in the diagonal rep- 
resentation is known one can compute the spectrum of 
the Hamiltonians. 

To find the eigenfunctions of the Hamiltonians, one has 
to solve a purely algebraic problem: one writes the one- 
dimensional representations in the basis of the regular 
representation. To find the ground state wavefunctions 
corresponding to the stationary states one has to con- 
sider the irrep in which all the generators are equal to 
one. The stationary states can easily be obtained for any 
Hamiltonian obtained from a DAA. The stationary states 
always have a product measure (there are no correlations) 
and depend on fewer parameters than the DAA. All the 
wavefunctions can probably be obtained. In a simple ex- 
ample (see Appendix A) we show not only how to obtain 
the whole spectrum but give also all the wavefunctions. 
They have a simple factorized form. 

We have studied some simple Hamiltonians in which 
only the generators appear. For a simple mathematical 
reason, in the finite-size scaling limit, they are all gapless 



with a dynamic critical exponent z — D (D is the number 
of dimensions of the lattice). There are many interest- 
ing time dependent physical properties which are worth 
studying in these models, the simplest being the time evo- 
lution of local densities. We plan to look at these aspects 
of the models in the future. Instead we have studied the 
properties of avalanches. Each DAA together with the 
ground state wavefunction of the Hamiltonian can define 
avalanches. 

In one dimension and in a DAA algebra which con- 
serves the number of particles, we have checked for two 
and three-state models that one gets avalanches in the 
" random walk" universality class pj . We have also stud- 
ied the "composition" of avalanches. One has a depletion 
of particles at the position where the avalanche starts and 
an enrichment at the end of the avalanche. The Hamil- 
tonian is gapless even if the DAA does not conserve the 
number of particles. This implies that one doesn't have 
SOC (the probability to have large avalanches decreases 
exponentially) in spite of the fact that the Hamiltonian is 
critical. One can complement this observation with an- 
other similar observation. Since we are in one dimension 
the stationary state of a A-state model can be seen as a 
spin model defined on the dual lattice. The spin model 
depends on many parameters (there is only one for N = 2 
corresponding to the temperature) and is not critical in 
spite of having a stochastic model which is critical. The 
opposite situation is already known (22J: a Hamiltonian 
having a gap drives the system to the critical (zero tem- 
perature) spin system. 

In two dimensions all we did is to check that for various 
choices of parameter one stays in the universality class 
described in [11, [II [11. 
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APPENDIX A: SPECTRUM AND 
EIGENFUNCTIONS OF A HAMILTONIAN 
DERIVED FROM A QUADRATIC DIRECTED 
ABELIAN ALGEBRA 

In section 2 and 3 we have discussed the quadratic 
abelian algebra 



eigenvalues of a^-i are: 



af = juaf +1 + (1- fi)aiOi +1 , i = 1,2,. 
a.L+i = 1, [a%, a,j] = 0, 

and the Hamiltonian: 

L 

L 



F = i]T(l- ai ). 



(Al) 



(A2) 



In order to obtain the spectrum of H we have to find 
the representations of Oj in the diagonal representation. 
In order to find the wavefunctions one has to express the 
dj in the regular representation. In the latter representa- 
tion, the vector space is given by monomials of the a,j 's 
(see section 2). 

In the representation in which a,i are diagonal, the dis- 
tinct eigenvalues Xi of a, can be obtained recursively from 
the eigenvalues x,+i of cn+\; 



{xi - x l+ i)(x,i + fix i+1 ) = 0. 
Equation (|A3[) can be written as: 

x, L = (- ( u) £i .T i+ i (a = o, l). 

This implies that the eigenvalues of ffl£_| are: 



,£L) 



IK-M) £ 

3=0 



(A3) 



(A4) 



(A5) 



and are labelled by the ej = Q, 1 (i = L, L — 1, . . . , L — 
I). In the 2 L representation of the algebra (|A1[) . the 



%L-l{Cl, £2) ■ ■ ■ , ZL-1-1,£L-1,£l) 
L-l-1 I 

= II 0-Y k tl(-ri ez -'> (A6) 

3=0 



k=l 



except for a±, they are degenerate. 

We are going to show that in the regular representa- 
tion, the factorized expression 



n 

m=0 



dL-r 



(A7) 



is an eigcnfunction of a^-i corresponding to the eigen- 
value (TA31): 



(A8) 



a L _ ; *(e L _ ; , . . . , ei_i, e L ) = x L -t(e L -i, ■ 
x*(e L _ ; , . . .,e L -i,e L ). 



This implies that \? eil ... )£L is an eigenfunction of all 
di (i = 1,...,L) corresponding to the eigenvalues (|A6|) . 
In particular if we take all the e, = 0, one obtains the 
eigenvalue one for all the a% which gives an eigenvalue zero 
for H (see (|A2I) ) and the corresponding eigenfunction is, 
up to a normalization factor, equal to (|16[) . 

The spectrum of H is obtained introducing the eigen- 
values (|A"6j) in (|A"2"1) . 

Obviously (IA7|) is also an eigenfunction for all a^-j 

(o < i < o- 

Using (|A6|) one can compute all the eigenvalues of H 
(they are fixed by the values of the e*, i = L,L — 1, . . . , 1). 
The corresponding eigenfunctions being given by (|A7|) 
with I = L — 1. The proof of (|A8[) is by finite induction. 

a) (1A"8)) is valid for I = 0. Using ([A"T|) one has: 

a L (Ai + (-/i) £i aL) = {-iiY L (}i + {-iiy L a L ), (A9) 



in agreement with (|A8[) . 

b) If (|AB|) is valid for i, it is valid for i + 1 We have 
to check that cll-i-i applied to ^(cl-i-i, ■ ■ ■) gives the 
eigenvalues (| A6|) . One has: 



= a i _ i _ 1 ( At +(-/i) £i - i - 1 ^=o 



~ 3 a L -i-i) 



x*(e L _i, 6l-j+i, ■ • • , e L ) = x L -i-x(ei,e2, 
x*(e L _;_i,eL_;, . . .,e L ), 



,£l) 
(A10) 



where we have used (|A1|) and (|A8|) . 

We think that the factorized form of the eigenfunctions 
(|A7[) is a general property of stochastic models obtained 
from DAA. It is a consequence of the representation the- 
ory of the DAA and not of a specific form of the Hamil- 
tonians. 
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APPENDIX B: SEVERAL REPRESENTATIONS 
OF THE QUADRATIC DAA ALGEBRA (5) 

We have to stress that unlike the spectra of evolution 
operators which are insensitive to the choice of the repre- 
sentations of the DAA algebras, the physical properties 
of the systems depend on the choice of the representation. 
For example, the existence or not of avalanches or even 
the universality classes of avalanches are not determined 
by the algebras but by their representations. Similarly a 
redefinition of the linearly independent words W(r) (see 
Eq. (JT)) ) in the algebra changes the physics. For example, 
the shift a 8 £ a- + 7 in the algebra (|10[) can make one 
of the coeficients P2 , Vz or 12 vanish, such a shift changes 
the physical process. 

In order to clarify the importance of choosing differ- 
ent representations for the same algebra, we consider the 
simpler case of the quadratic algebra (j25|) : 

af = fxa- +1 + (1 - fx)aiOi + i(i = 1,2,.. . L), 
aidj = a,ja,i, a L+1 = 1. (Bl) 

The expression of the generators in the 2 L regular rep- 
resentation used in Sections 2 and 3 can be written in 
the following way: 



with 



Ai 



(B2) 



where U,j is the 2°~ l unit matrix acting onto the tensor 
product Vi ® V(i + 1) ® ■ ■ ■ ® V(j), V is a 2-dimensional 
vector space. Ai is a 2 X 2 matrix, with elements acting 
on a 2 L ~ % vector space (order 1 and ai) 



Ai = 



nA 2 i+1 
h+i,L (1 - Li)A i+ i 



(B3) 



A second representation of the algebra is: 

a t = h ti -i <8> B t (B4) 

where 

' 1 - a (3 



V. = ( ~ a " 1 1 (3 ) ®B<+1, (B5) 



+ (3- 1. 



(B6) 



There are no avalanches in this representation. A particle 
hitting the site i doesn't influence the site i + 1. 

A third representation is used in the totally asymmet- 
ric Oslo model 101: 



ai = hi-! (8 Ci 



(B7) 



where 



V aI i+l,L (1 - 0)C i+ l J V ' 

with ji = a + /3 — 1. 

The eigenvalues of ai are the same for all the three 
representations. The eigenfunction <I> 



a 4 <I> = $ (i = 1,2, . . . ,L) 
is also the same: 



1 (P 



a + f3 \a 



(B9) 



(BIO) 



For the representation (|B3|) one has to take a = 1, /3 = 
/i in (|B9p and obtain (fTB)) which was used in Sections 2 
and 3. 

One can repeat the calculations done in Section 3 and 
derive the properties of avalanches. The results are a sim- 
ple generalization of the previous results. Equation |43|) 
stays unchanged with /i given by (|B6[) and the diffusion 
constant is D = a/3/ (a+/3) 2 . For a = 1, one recovers the 
results of Section 3. The parameter /i in the algebra (|Bip 
can take negative values. This is not the case if a = 1. 
The average density in the stationary state 



p = a/ (a + (3) 



(Bll) 



spans the interval < p < 1. The interval is smaller if 
a = 1. 



